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SUMMARY 


The effects which moderately Targe deflections have on the frequency 
spectrum of rotating rings and cylinders are considered. To develop the 
requisite solution, a variationally constrained version of the Lindstedt- 
Poincare procedure is employed. Based on the solution developed, in addi- 
tion to considering the effects of displacement induced nonlinearity, the 
role of Coriolis forces is also given special consideration. 

INTRODUCTION 

Numerous engineering applications (tires, turbines, satellites, etc.) 
contain rotor systems which are essentially rings or shells of revolution 
rotating about their axes. Obviously, in order to proper Ty influence their 
design^ a thorough dynamic analysis is necessary. In this regard, numerous 
papers have been published which deal with the free vibration properties 
of such systems. Most such work has centered on stationary configurations, 
as can be seen from the excellent surveys by references 1 and 2. The effects 
of rotation, in particular Coriolis forces, have been discussed by references 
3 to 7, With the exception of references 6 and 7 which treated small dynamic 
deformations superposed on large static deformations, the previous investiga- 
tions incorporating Coriolis acceleration forces have been limited to linear 
shell theories. This is a shortcoming since numerous rotor systems, tires, 
satellites, and turbines are flexible enough to undergo significant deflections 
in the form of moderately large rotations. 

It is the purpose of this paper to consider the effects which such 
moderately large rotations have on the frequency spectrum of rotating struc- 
tures. In particular, the analysis presented will consider the free 
vibration characteristics of rotating rings and cylinders wherein the deflec- 
tions involve moderately large rotations. Since the analytical model used 
to characterize the stated problem involves nonlinear partial differential 
equations, a modified version of the renormalized perturbation procedure is 
employed to evaluate the overall solution. This modification was undertaken 
since the usual renormalized procedure is unwieldy for systems of equations 
involving a multitude of frequency eigenvalue branches and secondly yields 
steady state results which are irregular for the linearized case. The modi- 
fication employed involves prescribing the system energy in advance; hence, 
a hierarchy of energy states is obtained from which the strained parameter 
can be evaluated. The resulting solution employing this procedure is 
regular, and thus, the proper limiting behavior is obtained for the linearized 
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case. Based on the solution, in addition to considering the global effects of 
nonlinearity, special emphasis is centered on determining the effects of 
Coriolis forces in the range of deformations marked by moderately large rota- 
tions. Hence the effects on the backward and forward traveling waves will be 
evaluated. 

GOVERNING EQUATIONS 

Since the nonlinear oscillations of rotating, elastically supported rings 
and infinite cylinders undergoing deflections involving moderately Targe rota- 
tions are considered herein, the governing displacement equations of motion 
employed to model the stated problem are defined by (refs. 2,4,6, and 7) 

|e2A^W,2w,Qg+f COS (me) cos (ut)+ph(W,^^-2fiVy^-n2W) = g (1) 


A^(V,ee+W,9+eW,0W,QQ)-ph(V,tt+2s^W,t-fi2V) = 0 

where 

A = — , A - — 

1 r4 2 r2 

such that e = W /R and 0, t, ( ),0, ( ),^, W, V, Wff,, E, I, h, R, p, P, k, w, and 
0 respectively ^ represent circumferential space, time, space and time dif- 
ferentiation, radial and circumferential shell displacements, maximum radial 
displacement. Young's modulus, moment of inertia, shell thickness, radius and 
density, internal pressure, foundation elasticity, exciting frequency, and 
lastly, the rotational speed of the shell. Due to the inherent nature of the 
circumferential coordinate space and the fact that the steady state response 
is being sought, it follows that W and V are periodic in both space and time. 

To round out the requisite field equations, the following potential energy 
functional is associated with equations (1) and (2), namely 

T 2 it 

Y = / / {A W,2^+A (V,2+2V,gW+W2)+eA {V,qW,2+WW,2) + 

0 0^2 2 

|e2A w,J+(ic+|)W2+2f cos (me) cos {a)t)W-ph[ft2(R*+W)^+W,2 + 

f22v2+V,2+2fi(R*+W)V,^-2j^W,^V]}dedt (4) 

where T = X— and R* = R/W„. 

CTTOi m 


( 2 ) 

(3) 
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SOLUTION 


As noted earlier, the standard renorraalized perturbation procedure has 
the twofold difficulty of yielding irregular results as e>0 and secondly, is 
unwieldy when more than one equation of motion involving several frequency 
eigenvalue branches is considered. This difficulty is circumvented by pre- 
scribing the systems potential energy in advance such that (W; V) = (W(8,t,f, 
m,Y); V(0,t,f,m,Y)). Once the solution is obtained, the role of y and W are 
reversed to that employed in the traditional version of the renormalized pro- 
cedure. To initiate the solution, o) is treated as the strained parameter; 
hence W, V, and w are expanded in the following perturbation series 

<W; V; o)> = I <W.; V.; a).>e’ (5) 

i=0 Til 

such that time is stretched so that x = wt. 


In order to obtain the zeroth order equations, e is set to zero; this 
yields 




-£^^W^)+f cos (me) COS (t) = 0 

(6) 

2 0,06 0,6' 0 0,TT 0 0,T 0' 

(7) 

2ir 2tr p 

r = / f {A W2 , +A (V^ +2V„ „W„+W2)+(k-*“)W2 + 

0 0 1 2 O>0 0,9 0 o' R' 0 

2f cos (me) cos (x)W -ph[fi2(R +W ) +0)2w2 +fi^V2+0)2V2 + 

0 0 0 OsX 0 0 o,x 

2a)„n(R*+W„)V^ - VJ}d0dx 

0 0 0,X 0 0,X 0 

(8) 

whereas with time, the potential energy space is stretched so that r = y/^. 
Since the steady state solution is sought. 

(“o’ ''o> ' <“co’ ''co> ''so> 

(9) 

where Wqo»... are time dependent. Employing equations (9), (6) 
to the following matrix set of ordinary differential equations, 

, and (7) reduce 
namely 

+[B,JY +f cos (x) = 0 
0 im .,mo,xx 0 2 m ^mo,x sm _mo ^ 

(10) 

such that 
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(11) 



(«co> ''so- “so- ''co>" 


Noting that [ 62 ^] is skew symmetric while [B^^] and [B ] are purely sym- 
metric, the steady state form of is given by . 


Y = Z cos (t) + Z„_ sin (t) 
^mo ^mc ^ms 


(12) 


where Z^^^ and Z^^^^ satisfy the matrix equation 



a)^[B,„]-[B,J 

0 im 3m 


oi T B ] 
0*- am-* 


0) fB 1 
0^ 2m-^ 


0) 




O'- im 



3m"'-J 'v^ms 


(13) 


Noting that the pencil of equation (13) yields the characteristic equation 
of equation ( 10 ), equation ( 12 ) becomes unbounded for coq equal to the natural 
frequency eigenvalues of the linear case. The properties of such eigenvalues 
can be ascertained by developing the appropriate Rayleigh quotient. This is 
possible by inserting' ’tmaZmeJT^ into^ equation ( 10 ) to yield a complex second 
order regular polynomial matrix problem. The inner product of this expression 
and Z^ yields a bilinear form from which the following modified version of 
Rayleigh's quotient is obtained, namely 


= J 



’&B3JZ^ I 


2 






V2 


(14) 


As can be seen from equation (14), Coriolis forces cause a twofold bifurcation 
in the number of eigenvalue branches. Following the previous comments, the 
relationship between r and wq, Wq, and Vg must be evaluated by inserting equa- 
tions (9) and (12) into equation ( 8 ); thTs yields 

[A^mHA^+K+l - ph(n^+io 2 m 2 )](w 2 ^^+w 23 ^+W 2 ^^+w 2 ^) + 

^'^sco'^'^cco ^ 

^sso'^cso''^cco*'^sco~'^cso*'^sso^"^^^“o^”’^*^cco'^cso"''^cso^cco ^ 


W V 
SCO sso 

'^sso^'^sco 
where W^^q,.. 


'''^sso'^sco''^cco'^cso^'^cso'^cco"''^sco*'^sso ^ 

)] = r/Tr2 (15) 

.,Vsso denote coefficients of the and solution, namely 
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^^ 0 ^ ^^cco {tti0)cos {t)+. . -+(^550 ;VggQ)sin (me )sin(T ) 

As can be seen from equations (13) and (15), four potential energy resonances 
are initiated for WQ~0(a)„,f) wherein are the frequency eigenvalues of the 
linear problem. Hence equation (5) is regular for e ^ 0 (the linear case). 

The first order set of field equations can be obtained by taking the first 
derivative of equations (1), (2), and (4) with respect to e and then setting e 
to zero. This yields 


^<k,e-''o.ee‘'o.6"''„,A,eo^“o«o.ee) ' 


2ir 2 tt 


(17) 


/ / 
0 0 




2(ic+^)WoWi+2f cos (ltl9) eos (T)Wj-2ph[!i2(R +W^|)WJ^ 
0 1 0,T 0 0,T 1,T 0101 0,T 

R V +Jl(o R V, W V 

0 0,T 1,T 1 0,T 0 1,T 0 1 0,T 

J2(i).W V V -Jlw.W V "ilca.W V - 

0 0 1,T 10 0,T 0 1,T 0 0 1,T 1 


fitOiWi^^Vo]}d0dT = 0 (18) 

Noting the form of the inhomogeneities appearing in equations (16) and (17), it 
follows that and V-j can be taken in the form 


(Wj;V^) = (W^,j;V^^) cos (2T)+(Wg^;Vg^) sin (2t) + 

•••■^^^ssi’'^ssi^ 

where the coefficients Wd,... are directly obtained upon inserting equation 
(19) into equations (16) and (17). Furthermore, employing equation (19) in 
conjunction with the first order potential energy constraint, equation (18), 
the following functional relationship is obtained for wi, namely 
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( 20 ) 


0 )^ = 0 )^ (1/D(2o)q, 0), 1/D(0,2m), l/D(2o)Q,2m) ) 

where D is the determinant of the pencil of equation (13). Hence for Wo~^(%f)» 
(01 is bounded and positive definite. This follows since l/D(2(Oo»0),. . .ete. 
remain bounded for V(Ooe(0,«>). Therefore, unlike the zeroth order set, Wi and 
Vi remain bounded for Vooo. 

In order to obtain the second order field equations, equations (1), (2), and 
(4) are differentiated twice and then e is set to zero. This operation yields 

= ''2[‘'u9«o.e"^.99“o.6■^''o.99^,9■^''l.9"o,99■^''o.e^.ee * 


«l‘'o,9"Vl.e9]4«.9“o,69-'>''[2“A«o,„" 

''2<''2,99"«2.9> - = 

■''2<“l.9“o,99'"'^0,9''l,96> * 

2»20Wo_+2»jSJWj_^) 


’>-i i <^«l9e*2^«„,99“2.e9■^^(^.9■^^)^ * 

2''2(''o.9-^Wo)(’'2,9"“2>*''2(2(''o.9"‘'o>'‘o.9Wl.e" 

2f cos (me) cos (t) - ph[ 2 j 22 (R*+WQ)W 2 +fi 2 w 2 + 

2W-(0 W ^ + 

O' 2,X 0,T 1,T 0 1 0,T 1,T 2 0 0,T 

^^(2V^V +V2)+u)2(2V, +V2 )+4(o u V, V„ + 

02 1 0 2,t 0,T 1,T 1 0 1,T 0,T 

2(1) (o,V2 +2J2(»).(R +W )V +2fiu) (R +W )V + 

0 2 o,T 0 0 2,T l' 0 1,T 

2J2(o (R +W )V +2s^o) W V +2^(1) W V +2J2 co_W V - 

2' O' 0,T 0 1 1,T 1 1 0,T 0 2,t 0 


( 21 ) 


( 22 ) 


(continued) 
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(23) 


2 fio) W V - 2 no) W V.-2flco W V - 2 to W V 

0 0,T 2 0 2,T 0 1 0,T 1 1 1,T 0 


2 nw W V ]}d6dT 
2 0,T 0 


As in the zeroth and first order cases, noting the inhomogeneities of 
equations (21) and (22), W2 and V2 take the form, namely 

cos {Zt) + 


+ (Wss2’ '^SS2^ cos ( 3 t) ( 24 ) 

Employing equations ( 24 ), ( 21 ), and ( 22 ), it can be shown that the following 
proportionalities exist, that is 


(W^; V2)“(W2(l/D3(wQ,m)); 


( 25 ) 


Hence W2 and V2 become unbounded for a3p~0(a)mf). The requisite form of C02 can be 
obtained by inserting equation ( 24 ) into the second order potential energy 
functional, namely equation ( 23 ), After extensive manipulations, this operation 
yields the following proportionalities for U2, that is 


'^2 ^ 


( 26 ) 


Thus for too~ 0 (o)n)f) , a)2~0(l/D2(a)o>m)) where, since D^(wpif,m) is singular, 0)2 
is itself unbounded and negative definite. Additionally W2 and V2 are 
themselves unbounded at such values of coq. 


DISCUSSION 


Stopping the solution at this point, W, V, and w are given by 
(W; V; oj)~(Wq; V^; w^) + (W^; V^; (o^) e + 

(W^; y^i + 0(e3) (27) 

Due to the procedure employed, it follows that W and V are regular in e, 
including s = 0 . This result is, in contrast to standard renormalized 
perturbation procedures which do not yield zeroth order solutions exhibit- 
ing the proper unbounded behavior for u on the order of the linear system 
frequencies. 

The softening behavior of the ring or infinite cylinder can be directly 
obtained by considering the fundamental relationship between 03 and r. 

Before doing this, the nature of the wq dependency of w must be ascertained. 
In particular, for a)o~ 0 (a)nif). 
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u)~03 + eO(1) - e^O (■ — -L- — )+ O(e^) (28) 

° 02(0)^, m) 

where since (i >2 is negative definite and unbounded, m is itself negative 
definite and unbounded. Such unboundedness occurs at each of the eigen- 
values of the pencil of equation (13), Note as Q is set to zero, the tw6 pairs 
of eigenvalue branches merge back to the two frequency branches of the station- 
ary state, and hence, the traditional frequencies are obtained. 

Eliminating cog from equations (28) and (15), it follows that since w is 
unbounded and negative definite for «o-0(cqpf), the overall steady state harmonic 
behavior of the ring or infinite cylinder is of the softening type. Hence, 
as w is raised or lowered, the usual softening type jump phenomenon is 
encountered. 

In the context of the foregoing, the results can be summarized by the 
following remarks: 

(1) Coriolis forces induce bifurcations in the frequency spectrum; 

(2) Such bifurcations extend into the range of deflections marked by 
moderately large rotations; 

(3) All branches exhibit a softening type behavior; this applies to the 
branches associated with forward as well as backward traveling waves; 

(4) Driving frequencies in the neighborhood of the linear system 
frequency may induce jump phenomena; 

(5) Setting fi->0 yields the results for stationary rings and cylinders. 
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